Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Massive zero-metal stars: Energy production and mixing 
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Abstract. Time-dependent nuclear network calculations at constant temperature show that for zero-metal stars 
> 20Mq (i) /9-decay reactions and (ii) the ^^N(p,7)^^0 reaction must be included. It is also shown that the 
nuclear timescale in these zero-metal stars is shorter than the mixing timescale and therefore the assumption 
of instantaneous mixing across convective regions is not fulfilled. We conclude that proper modeling of these 
processes may alter the evolution of massive zero-metal stars. 
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1. Introduction 

The evolution of zero-metal stars in the mass range ~ 20- 
IOOMq , i.e., massive Population III stars (Pop-Ill for 
short), has been studied since the pioneering work of Ezer 
|et al. (1971 ) who followed the evolution of stars in the 
mass range 5-100Mq from the pre-main sequence con- 
traction phase until the exhaustion of hydrogen on the 
main sequence. This opened the field for many studies 
t hat were mainly concerned w i th later evolutiona r y stages 
(ICary, 1974 [Castellani, 1983| ; |E1 Eid et al., 1983t jober ei " 
|al., 1983| ; |Klapp, 1983| , |1984| ). Motivated by the still ongo 



ing debate about the initial mass function of Pop-III stars, 
pVIarigo et al. (2001 ) present the most comprehensive study 
of zero-metal evolutionary models (0.7-100Mq) starting 
from the ZAMS until the AGB in the case of low- and 
intermediate-mass stars, or to the onset of carbon burn- 
ing in massive stars. 

In earlier studies the authors were forced to make nu- 
merous assumptions about equilibria between the chemi- 



cal species involved. The most recent studies by Marigo el 



al. (2001), however, incorporate nuclear networks without 
making any assumptions about equilibria. Nevertheless, 
even these recent state-of-the-art studies rely on some sim- 
plifications that are certainly fulfilled during the evolution 
of normal stars but are questionable in the case of massive 
zero-metal stars. 

It is the aim of this letter to show that the following 
two assumptions widely used, namely 

1. /3-decay negligible against proton-capture 

2. instantaneous mixing across convective regions 



are not met during the evolution of massive Pop-III stars 
and that proper modeling may alter the evolution of these 
objects. 

2. Nuclear Network 



It was first noted by Ezer (1961) that Pop-III stars more 
massive than about 20Mq are supplied by a different 
mode of energy generation. Due to the lack of the elements 
C, N and O the CNO-cycle cannot generate the luminosity 
needed to halt contraction. This raises the temperature to 
the point where the 3a-reaction produces enough carbon 
to initiate the CNO-cycle. As a result, the CNO-cycle op- 
erates at a much higher temperature, typically at 10^ K. 
Since the energy generation in the CNO-cycle exhibits a 
strong temperature dependence the central temperatures 
of the more massive mo dels up to 100M(^ do not exceed 



Tc = 1.3 10^ K 

[2001 ). At these temperatures 



(see e.g. 



El Eid et al., 1983; Marigo et al 



the CNO-cycle turn-over 



rate is governed by the beta-decay half- lives ( Ezer et al. 
19711) . 



This finding has one important consequence for mod- 
eling stellar evolution: one has to include the elements "'^^N 
and -^^O (mono-cycle), ^^F (bi-cycle) and (tri-cycle) 
explicitely in time-dependent network calculations. In ad- 
dition, for temperatures exceeding 10^ K the proton cap- 
ture i3n(p,7)14o (followed by "0(/3+,i/)"N) supersedes 



theJ£N(/3+,i/)^^C rate. Contrary to the findings of Klapp 
1983D our calculations show a non-negligible effect on the 
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energy generation for T > 10* K. 

In order to estimate the described effects, we have per- 
formed time-dependent network calculations at constant 
temperature. Our reference calculation includes 15 chem- 
ical species: ^H, ^He, ^^c, ^^C, ^^N, "n, i^n, i^q, isq. 
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Fig. 1. Energy generation per pX^ at constant temperature for a) T = 1.0 lO^K and b) T = 1.3 lO^K. The solid curve 
is our reference calculation including all 15 chemical species. The dashed curve neglects the /3-decay and the dotted- 
dashed curve neglects the reaction ^'^N(p,7)^^0. Differences between the reference curve and the curve neglecting 
^'^N(p,7)^''0 are too small to show up in a). 



1^0, "O, 180, "F, 
of Xh = 0.77, Xne 



i^F and ^^F with initial abundances 
= 0.23 and Mother = 10^1^, and the 
following reaction chains: 



- 3Q;-process 

- CNO tri-cycle 

- i3n(p,7)140 and ^^0{(3+ ,vf^n 



Reaction rates are taken from A^j4Ci?£'[| (see Angulo et al 
1999|)~|The calculations are performe d utilizing the DA E 
solver LIMEX that is maintained by Ehrig et al. (199^ ) . 

Our results can be viewed in Fig. |l] where we have 
plotted the energy generation per pX'^ against the time in 
years for the two cases T = 1.0 10^ K (see Fig. |l|a) and T — 
1.3 10^ K (see Fig. ^). Neglecting the beta-decay half-lives 
(dashed curves) has an impact on the energy generation 
of ~ 10% for 1.0 10^ K and - 65% for 1.3 10* K. Whereas 
the proton capture on I'^N alters the energy generation at 
1.0 10® K only slightly (< 1%), it has a - 30% effect at 
1.3 10® K (dotted-dashed curve). It is worth noting that 
the energy generation is only marginally (< 1%) effected 
by both the second (named ^^ON) and the third cycle 
(named ^^FO). These cycles need only be included if one 
is interested in the relative abundances of the elements O 
and F. At constant temperature, the energy production 
never attains a constant value due to the ongoing feeding 
of freshly produced carbon via the 3a-process. 

3. Mixing 

Massive stars contain large convective cores in which mix- 
ing of chemical species occurs due to the turbulent con- 
vective motion. In normal stars, this process is very rapid 
compared to the slow changes of the chemical composition 
produced by nuclear reactions. Under these circumstances 
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one can safely assume that the composition in a convec- 
tive region always remains homogeneous, i.e., elements are 
instantaneously mixed over the convective region. The fol- 
lowing arguments are put forward to show that the turn- 
over time of the CNO-cycle at the high temperatures of 
massive Pop-Ill stars can be comparable to the timescale 
of turbulent mixing - or even shorter. 

3.1. Mixing timescale 

We show that the mixing timescale, Tconv, is of the or- 
der of 10 days in the case of both normal and zero-metal 
stars. An estimate of this timescale is easily derived from 



the standard mixing length theory (MLT, Bohm-Vitense 



1958) and the formulas we present resemble those given 



in Kippcnhahn (1990, chap. 7). Let us assume a standard 
composition of = 0.77 and Xhc = 0.23, i.e., — 0.58, 
and opacity due to electron scattering k — 0.2(1 + Xh), 
i.e., n = 0.35cm^g^i. Even in the case of hot Pop-Ill 
stars, neglecting radiation pressure does not change this 
order of magnitude estimate and it suffices to assume a 
monoatomic ideal gas: 5 = 1, cp = 51K/2/i. The mix- 
ing timescale is approximately given by the pressure scale 
height 77p divided by the velocity of convective motion 



With this: 



10** 



M, 







T f IQ- 



(1) 



^RqJ \ m J \I 10®K VVV - Vc 

where r, T and m are the radius, temperature and mass, 
respectively, at one particular locus in a star. The quantity 
V — Vo is the difference between the actual temperature 
gradient, V, and Vc which describes the variation of T in 
a mass element during its motion. Introducing two dimen- 
sionless variables: 

3/2 



[/ ~ 210" 
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and provided that ?7 < (V - Vo)^/^ < W (see below) the 
cubic equation of MLT can be simplified yielding 



1/3 



(4) 



For an estimate of Tconv at the center of a star the first 
two equations can be even further reduced, since m = 
4/37rper3: 
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(5) 



and 

f7 ~ 410' 



10 



M, 
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(6) 



For a typical example, consider a 20Mq Pop-III star. 
There we have Tc ~ lO^K, pc ^ 100 gcm~'^. Let us choose 
a mass coordinate: to = I-Mq, hence U ~ 410"^" from 
Eq. d). For reasonable values of M^, - 1 ... 100 it fol- 
lows that (V- Vo)^/^ 7 10^^. Note that the assumption 
U < (V- Vc)^/2 < is fulfilled. With Eq. (|) we finally 
arrive at 



10 days . 



(7) 



For a star containing metals the central temperature is 
by a factor of three smaller but as normal stars are less 
compact also the density goes down by a factor of ten. 
Thus Tconv changes by a factor of two or so. Similar results 
are obtained at other loci of the convective core and for 
stellar masses in the range 20-100 Mq. 

3.2. Nuclear timescale 

A characteristic timescale for nuclear reactions in massive 
stars is the turn-over time of the CNO-cycle. This time is 
primarily given by the slowest reaction of the cycle which 
is ^''N(p,7)^'^0. Even in the high temperature regime (T = 
1.3 lO^K) where proton-captures become comparable to /3- 
decays this estimate holds within a factor of two. Hence, 

Ah 



'nl4pg 



-^HpAnl4pg(7') ' 



(8) 



where Ah, and p are the atomic mass of hydrogen, 
the mass abundance of hydrogen and the density, respec- 
tively (see playton, 1983 ). The reaction rate Ani4pg(r) 
measured in [cm^^ Mol^-'- s^-'-] is only a function of temper- 
ature. Taking this reaction rate again from NACRE (see 



Angulo et al., 1999) 



A„i4pg(T) = 4.83 10^ Tg 



-2/3 
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Fig. 2. Lines of constant Tnuc at Xh — 0.77. In the region 
above the dashed line the mixing time is slow against the 
nuclear timescale. 



X exp -15.231 Tg - ( 77^ 




X (1 - 2.00 Tg + 3.41 - 2.43 T^) 



2.36 lOTc 



-3/2 



g exp(-3.010rg"') 
6.72 lO^TgO-sso exp(-9.530rg~i) , 



(9) 



the nuclear timescale is easily calculated for a given set of 
density p and temperature T. To explore the parameter 
space we plot lines of constant time Tnuc in a p— T diagram 
(see Fig. 2). The dashed line marks the typical timescale 
of convective mixing. Thus, the region above (below) this 
line is characterized by slow (fast) mixing compared to 
the nuclear timescale. The assumption of instantaneous 
mixing is only justified in the lower region where mixing 
is fast. It can be nicely seen that all normal stars lie in the 
lower region but that massive Pop-III stars lie in the upper 
region where the assumption of instantaneous mixing does 
not hold. 



4. Conclusions 

Performing time-dependent nuclear network calculations 
at constant temperature we show that in the case of 
massive Pop-III stars (> 20 Mq) neglecting the /?-decay 
against proton-capture leads to a considerable error in 
the energy generation rate. In addition, the reaction 
^^N(p,7)-'^^0 cannot be omitted for temperatures exceed- 
ing 10^ K. 

Moreover, the nuclear timescale of massive Pop-III 
stars can be very short (order of hours) compared to the 
timescale of convective mixing which is of the order of 10 
days. Therefore, instantaneous mixing which is well justi- 
fied in normal stars may introduce large errors in evolu- 
tionary calculations of zero-metal stars. 

Contrary to normal stars massive Pop-III stars have 
not forgotten their nuclear history since equilibrium 
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abundances between chemical species are never attained. 
Hence, evolutionary calculation must start on the pre- 
main sequence well before the onset of nuclear reactions. 
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